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Solutions to the complementarity problem constructed in [1], generally, pos- 
sess non-zero total charge. In natural sciences, bodies possessing non-zero to- 
tal charge (ions and similar object) are considered as specific objects. Bodies 
possessing zero total charge (electro-neutral bodies) are considered as general 
case objects. This paper presents a solution to the complementarity problem 
for electro-neutral bodies. The solution is constructed under the condition 
that the volumes of the bodies are small. 



The statement of the complementarity problem for charged bodies 

Let us consider two bodies with axisymmetric distributions of charges 0(x) 
and ^>(x) as shown in Fig.l. In this case, the interaction of two bodies 
is measured by the interaction force in the direction Ox^ [1]. We restrict 
the translational degrees of freedom in Ox\- and Oa^-directions and allow a 
translation along the Oa^-axis only 




Figure 1: A pair of bodies with axisymmetric distributions of charges 
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The total interaction force / (0, ip) for the given pair of bodies with dis- 
tributions of charges 0(x) and over them in the "one touches another" 
position is [1] 

1(0,^)=/ / J R(|x-y|)0(x)dx^(y)dy (1) 
Let us introduce the operator [1] 

R : 0(y) G L 2 (Q) -+ [ flflx - y|)0(y)dy G L 2 (Q). (2) 

JQ 

With this notation, we have the following formula for the pair interaction 
force 

/<0,V> = (R0,V), ( 3 ) 

where (•, •) means the standard scalar product in L 2 (Q) [2]. 

The kernel i?(|x — y|) of the integral operator in (2) is the fundamental 
solution to the electrostatic problem and is singular [3]. We make a simpli- 
fication and assume that the kernel is a smooth function (thus we assume 
that the charges belonging to one body do not touch the charges belonging 
to another body). 

A body is called electro-neutral if its total charge equals zero. In other 
words, the distribution 0(x) has zero average value over the domain Q 

(0) = / 0( x )dx = 
JQ 

(here one should not mix the averaging operator (•) with the coupling oper- 
ation ( , }). 

The subspace L = {0 G L 2 (Q) : (4>,ip) = 0} is the space of charges of 
electro-neutral bodies. 

Generally speaking, the operator R (2) acts from L to L 2 (Q), but not to 
L [1]. In other words, operator R transforms an electro-neurtal distribution 
into a distribution, that is not necessarily electo-neutral. 

Definition. We say the distributions of charges 0(x) and ip(x) to be com- 
plementary ifl{<f>,ip) > (4>(x) andifj(x) attract one another) and to be not 
complementary if I (<p,ip) > (<f>(x.) and "0( x ) push away one another). 

The problem is formulated as follows: "Do there exist four distributions 
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1> 





which belong to Lq and possess the following properties: is complementary 
to $ and not complementary to 0, ip and ty; is complementary to \1/ and 
not complementary to 0, ip and $, neither $ is complementary to \&?" 

This means that attracts $ and pushes away the others and vice versa, 
ip attracts \& and pushes away the others and vice versa. In all other possible 
pairs, the distributions push away one another. It is seen that the com- 
plementarity problem is a version of the "lock and key" problem [4] under 
the additional condition that the lock "attracts" a suitable key and "pushes 
away" all the others. Due to (3), the conditions above may be written in the 
form 

J(0,$)>O, (4) 
7(0,0) <O,/(0,V) < 0,1(0,*) <0, !<<£>,<£>) <0, <0, 

I >0, 

i {if;, ij) < o, / (v>, 0) < o, / (v>, $) < o, i *> < o, / $) < o. 

Due to the symmetry of the interaction forces with respect to the distribution 
of charges (see formula (1)), only ten inequalities in (4) are independent. The 
system of inequalities (4) must be solved under the condition 0, -0, $, * e L . 

Solution to the problem 

We have already noticed that L is not a proper subspace of the operator R 
(2). As a result, we cannot apply the method developed in [1] to solve the 
complementarity problem for electro-neutral bodies. 

We note that L is a subspace of L 2 (Q). Thus, L 2 (Q) can be represented 
as the sum: L 2 (Q) = Lq®Lq of L and its orthogonal complement Lq [2]. In 
the case under consideration, Lq = {0 e L 2 (Q) : 0(x) = const}. We denote 
by Pr the operator of projection [2] acting from L 2 (Q) to L . Then 

Pr0 = 0(x)-(0). (5) 

The linear operator R : L 2 (Q) — > L 2 (Q), which is compact and self- 
adjoint [1], has a system of eigenfunctions {(Pi}^, which form an orthogonal 
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basis in L 2 (Q). In particular, 



(6) 



where Aj means the eigenvalue corresponding to the eigenfunction ipi. 

We consider two functions Prfa = fa — (fa) and Prfa = fa — (fa) and 
compute the interaction force for this pair 



I(Rfa,Rfa) = (RPr^Pr^-)- 
By using (5) and (6), we have 

(RPiY?j, Pr<£j) = 

(R(fa-(fa)),fa-(fa)) = 
= (Xifa - (fa) Rl, fa - (fa)) = 
= Kifa, fa) - (fa) (fa, Rl)) - \i (fa) (fa) + (fa) (Rl) (fa) = 
HSij - (4>i) (fa)) - (0.) «0„R1) - (Rl) (fa))- 

We use here the equality (/, C) = C ■ (/), where C is a constant. 



(7) 



(8) 



With regard to the Schwartz inequality 



/(x)c/(x)rfx 



Q 



< 



Ml [2], 



we have the following estimates for the terms forming the r.h.p. of (8): 



0i(x)l<ix 



a/ mesQ, 



(9) 



(fa,Rl)\<\\fa 



IRlll = IIRll 



In (9), we use the fact that \\fa\ \ = 1 (where || || = 1 means the standard 
norm in L 2 (Q) [2]). Then the terms in the r.h.p. of (8) can be estimated as 
follows: 



lA^y - (fa) (fa)) - (fa) ((faRl) - (Rl) (fa))\ < 



(10) 



< mesQ + -\/mes(5(||Rl]| + |R1| a/ mesQ). 



By definition, 



Rl 



i2(|x-y|)dy. 



(11) 



We consider a set of domains Q(r) depending on the real parameter r and 
embedded one into another in such a way that Q(ri) D Q(r 2 ) if r\ > r 2 . Then 

Rl(r) = / i?(|x — y|)dy, being considered as a function of r, decreases. 

jQ(r) 

Then Rl(r) < / i?(|x — y\)dy = C for any r < r and the r.h.p. of (10) 
is not greater than 

mes Q + \pmes~Q • Cq{\J mes Q + mes Q^/mesQ) < C ■ mes Q (12) 

for any r < r , where C < oo. Thus, we can write (8) as 

(RPr^, Pr^) = XiSij + F tj , (13) 

where \Fij\ < C ■ mesQ, see (12). In order to construct four complementary 
functions, we consider the functions displayed in Table 1. 



Table 1. Two pairs of electro-neutral distributions 



<p = Pry?; + aPnp k 


$ = —Pripi + aPnp k 


ip = Pr(pj + aPrip k 


\I/ = — Pripj + aPrifk 



where ipi, ipj and ip k are eigenf unctions of the operator R, i ^ j, i ^ k, j ^ k. 
The functions </>, <£>, ip and \l/ determined by Table 1 belong to L (all these 
distributions are electro-neutral, since all of them are sums of projections of 
the corresponding functions to L ). The interactions corresponding to the 
functions displayed in Table 1 are the following: 

I. The interaction of <p with $ 

(R<j>, $) = -(R(Pr^ + aPnp k ), -Pr^ + aPr^) = (14) 
= (RPr<^, Pr^) - a(RPnp k , Pry^) + 
+a(RPr(pi, Pnp k ) + a 2 (RPry? fc , Pry?*.) = 
= -Xi - Fa - a(F ki - F ik ) + a 2 (X k + F kk ) = 
= -Xi + a 2 X k + [-Fa - a(F ki - F ik ) + a 2 F kk \. 

II. The interaction of ip with \I> 

(R^, *) = (R(Pr^ + aPry? fc ) , -Pr^- + aPr Vk ) = (15) 
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= -(RPr<£j,Pr</?j) - a(RPiY> fc ,Pr<£>j) + 
+a(RPr^, Prv9 fc ) + a 2 (RPry9 fe , Pr^ fe ) = 
= -Xj - Fjj - a(F kj - F jk ) + a 2 {X k + F kk ) = 
= -Xj + a 2 X k + [-Fjj - a{F ki - F ik ) + a 2 F kk \. 

III. The interaction of (f) with ip 

(R0, V) = (R(Pr^ + aPnp k ), Pnpj + aPnp k ) = 
= (RPr^j, Pripj) + a(RPry? fc , Prpj) + 
+a(RPrc^, Prip k ) + a 2 (RPr^ fc , Prip k ) = 
= Fa - a(F kj + F lk ) + a 2 (X k + F kk ) = 
= a 2 X k + [Fij - a(F kj + F ik ) + a 2 F kk }. 

IV. The interaction of with \1> 

(R<£, *) = (R(Pr^ + aPr<p k ), -Pr<^ + aPnp k ) 
= (KPripi, Pr<pj) - a(RPr^, Pr^-) + 
+a(RPripi, Pr(p k ) + a 2 (RPry? fc , Prtp k ) = 
= F i3 - a(F kj - F ik ) + a 2 {X k + F kk ) = 
= a 2 X k + [F^ - a(F kj - F lk ) + a 2 F kk }. 

V. The interaction of $ with ^ 

(R$, *) = (R(Pr^ + aPry? fc ), -Pr<pj + aPr^) 
= -(RPr^j, Pnpj) - a(RPrv? fc , Pr<^) + 
+a(RPrip h Pr(p k ) + a 2 (RPry? fc , Pnp k ) = 
= -F^ - a(F kj - F ik ) + a 2 {X k + F kk ) = 
= a 2 X k + [Fij - a(F kj - F ik ) + a 2 F kk \. 

VI. The interaction of ip with <£> 

(R^, $) = (R(Pr^ + aPry9 fc ), -Pr^ + «Pr^) 
= -(RPryjj, Pry?;) - a(RPrv? fc , Pr<^) + 
+a(RPr^-, Pry? fc ) + a 2 (RPry? fc , Pr^ fc ) = 
= Fji - a(F kl - F jk ) + a 2 (X k + F kk ) = 
= a 2 X k + [F^ - a(F ki - F jk ) + a 2 F kk }. 
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VII. The interaction of with </> 



(R0, 0) = (R(Pr^ + aPr<p k ) , Pr<ft + aPnp k ) = (20) 
= -(RPry?;, Pry?;) - a(RPry9 fe , Pry?;) + 
+a(RPry?j, Pry? fe ) + a 2 (RPrip k , Pry? fc ) = 
= X t + Fa - a(F kt - F lk ) + a 2 (X k + F kk ) = 
= X t + a 2 X k + [Fa - a{F kl - F lk ) + a 2 F kk }. 

VIII. The interaction of with ■?/> 

(R0, V) = (R(Pr^- + aPnp k ) , Pr^ + «Pr^) = (21) 
= — (UPrtfj, Pr<pj) — a(RPr(p k ,Pr(fj) + 
+a(RPr^, Pry? fc ) + a 2 (RPry? fe , Pnp k ) = 
= Xj + F j3 - a(F kj + F jk ) + a 2 (X k + F kk ) = 
= Xj + a 2 X k + [Fjj - a(F kj + F jk ) + a 2 F kk \. 

IX. The interaction of $ with $ 

(R$, $) = (R(-Pr^ + aPry? fe ) , -Pr^ + aPi%) = (22) 
= (RPry?;, Pry?;) - a(RPry? fc , Pry?;) - 
-a(RPry?;, Pry? fc ) + a 2 (RPry? fc , Pry? fc ) = 

= Xi + Fa - a(F ki + F ik ) + a 2 (X k + F kk ) = 
= Xi + a 2 X k + [Fa - a{F kl + F lk ) + a 2 F kk }. 

X. The interaction of \1/ with ^ 

(Rtf , *) = (R(-Pr^ + aPry? fe ), -Pr^ + «Pr^) = (23) 
= (RPry? j5 Pry?j) - a(RPr^, Pry?.,-) - 
-a(RPr^j, Pry? fc ) + a 2 (RPry? fc , Pry? fc ) = 
= X 3 + Fjj - a(F kj + F jk ) + a 2 (A fc + F kk ) = 
= Xi + a 2 X k + [F^ - a{F k j + F jk ) + a 2 F kk ]. 

Since the operator R is negatively determined [1], all of its eigenvalues 
{Xij^Li are negative. 

The r.h.p.'s of formulas (14)-(23) involve the following independent pa- 
rameters: the real number a and certain combinations of the quantities {i 7 ^} 
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with various indices (see the expressions in square brackets). The absolute 
value of the terms in the square brackets in (14)-(18) is not greater than 



F + 2aF + a 2 F < F(l + afmes Q, 



(24) 



where F = max{F mn ;m,n = i,j,k}. For a < 1, the r.h.p. of (24) is not 
greater than AFmes Q. If we take 



then the r.h.p. 's of (14) and (15) (the interactions of <ft with <3> and of ip 
with \l/) are positive and the r.h.p. 's of formulas (16)-(23) (i.e. all the other 
interactions of four distributions of charges displayed in Table 1) are negative, 
thus (4) is satisfied. 

This means that there exist two pairs of electro-neutral distributions of 
charges: one pair is cf) and $, another pair is ip an d \P and all the other 
possible pairs of distributions are not complementary. 

Inequality (26) means that the volume of the bodies is bounded from 
above. Note that although the volume mes Q is small, it is a finite positive 
number. 
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